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\Q ■ A model calculation is given for the energy relaxation of a non-equilibrium distribution of hot 

electrons prepared in a metallic sample that has been subjected to homogeneous photo-excitation 
by a femtosecond laser pulse. The model assumes that the delta pulse photoexcitation creates two 
CNj , interpenetrating electronic subsystems, initially comprising a dilute energy-wise higher-lying non- 

• degenerate hot electron subsystem, and a relatively dense, lower-lying electron subsystem which 

' is degenerate. In the femtosecond time regime the relaxation process is taken to be dominated 

by the electron- (multi) phonon interaction, resulting in a quasi-continuous electron energy loss to 
the phonon bath. The kinetic model is given for this time regime, beacuse in this time regime 
the usual Two Temperature model is not applicable. The Two Temperature model assumes that 
the hot electrons and phonons are in their respective equilibrium states (Fermi and Bose) but at 
a different temperatures. One uses the Bloch-Boltzmann-Peierls transport formula to calculate the 
energy transfer rate. In the present Kinetic model a novel physical feature of slowing down (due 
to Fermionic blocking of interaction phase space) of electron-phonon relaxation mechanism near 
. the Fermi energy of degenerate electronic subsystem is considered. This leads to a peaking of the 

calculated hot electron distribution at the Fermi energy. This feature, as well as the entire evolution 
of the hot electron distribution, may be time-resolved by a femto-second pump-probe study. 
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PACS numbers: 71.10.Ca,71.20.Nr,72.15.Lh 



I. INTRODUCTION 

s 

' Consider the kinetic evolution of a photoexcited nonequilibrium electron-phonon system, comprising three subsys- 
tems. The first two are interpenetrating electronic subsystems, one degenerate and the other nondegenerate, while 
the third is the phonon bath (the lattice). For these three interacting subsystems, one can define three characteristic 
O \ relaxation times, namely, the electron-phonon (r ep ), the phonon-phonon (t pp ), and the electron-electron (r ee ) relax- 
ation times. Thus, for example, in the case of metals, one usually has r ee <C r ep 3> r ep . In fact, the electron-electron 
' relaxation time is approximately one-tenth of that of the electron-phonon, the latter being a few picoseconds. For 
^ . this very reason, one uses the standard Two Temperature model^^i^. The present model can be used for the initial 
femto-second time scales where the electron distribution is still in non-equilibrium (non Fermi-Dirac Distribution) 
■^J" , state. 
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Consider a delta pulse photoexcitation creates two interpenetrating electronic subsystems, initially comprising a 
dilute energy-wise higher-lying non-degenerate hot electron subsystem, and a relatively dense, lower-lying electron 
subsystem which is degenerate. In the femtosecond time regime the relaxation process is taken to be dominated by the 
electron-(multi)phonon interaction, resulting in a quasi-continuous electron energy loss to the phonon bath. The fem- 
tosecond photoexcited non-degenerate hot electron distribution then evolves predominately through the non-radiative 
electron-phonon processes. In the present work, two models of electron- (multi) phonon interaction are considered. A 
linear model in which the phonon friction is taken to be linear in velocity of the hot electrons, and the other, in which 
the phonon friction is nonlinear in electron velocity. The nonlinearity considered here is a power-law type. In the 
following, we drive analytical expressions for the time-dependent non-equilibrium hot electron distribution for the two 
models of electron-(multi)phonon interaction. 



II. THE MODEL 



A. The phonon friction is proportional to the electron velocity V: linear model 



The kinetic equation for electron energy loss is: 

df e (e,t) . df e (e,t) _ df e (e,t) 

&t de r r " { ' 

Here,e = (mv 2 )/2 = —"f p v 2 — — e/r p , where r r and r p are the radiative(photon) and the non-radiative (phonon) 
relaxation times, and r = t/rp and a r _ p = T r /r p are dimensionless variables/parameters. With this, equation (1) 
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FIG. 1: Pile-up of hot electrons near the Fermi energy (r? = 0.5). 



reduces to, 

9/ e (e,r) df e (e,T) df e (e,r) 



(2) 



dr dine a r - p 

Equation (2) can be solved analytically to give: 

/ e (e,r) =f e- T '°"->e(-T-ln[e/(e F + e L )]) (3) 

Here ©[...] is a step function, and we have assumed an initial delta- function laser pulse of photon energy ej, that 
excites electrons in the energy interval ep < e < ep + cl- Inasmuch as Eq.(2) has only the forward propagating 
solutions (in energy space) we can readily incorporate the slow process effective at the bottom of the conduction band 
by simply introducing a longer relaxation time T pe > t p there. This gives rise to a kind of peaking effect, beacause of 
a pile-up of the hot electrons just above the Fermi energy of the degenerate electron subsystem. This peaking is due 
to the fermionic blocking of the electronic transitions, and it restricts the interaction phase space. Equation (3) is 
our basic result. In terms of it, we can calculate the total number N pe (r) of the hot electrons piled up near the Fermi 
energy, and also the total number of hot electrons A r ^ ot (r) above ep. The number of hot electrons in the pile- up, 
Af>(r) = N hot (r) — N pe (r), with energy e > e F (for uniform excitation of scale height, / ) is: 



dN pe (r) cW>(r) dN pe {r) 



/V = ^ (4) 



Now ^V>(t) = J^ +£L f e (e,T)de — foe T / a '—p J^ eL Q(—t — ln(e/[ep + ep])de. We consider energy loss through 
phononic friction only, with this we get 

JV>(r) =h{e c + e L )e- T (5) 

From equations (3) and (4) we have 

- fo(e F + e L )e^ = - V N pe (r) , V = ± (6) 
or p pe 



Solving by Laplace transforms we have 



N pe (r) = / ° (£F 1 + _^ )e " T [1 - e-^l (7) 



and the total number of hot electrons 



N hot (r) - N pe (r) + 7V>(r) = f (e F + e L )e~ T 



1 + - (1-e-^) 

1—7/ 



(8) 



Time evolution of these two populations are plotted in Figs.(l) and (2). 



3 



Nhot ( -O ) 




7T 



2 4 6 8 1 O 

FIG. 2: Decay of hot electrons, N hot (T), for 77 = 0.5. 



B. The phonon friction is nonlinear (algebraic) in the electron velocity: nonlinear model 



dt 

With e = — ^[^]™ and £o is an associated energy scale. 



df e (e,t;n) df e {e,t;n) _ f e (e,t;n) 



de 



(9) 



In line with the linear model, this kinetic equation has the solution, 



/ e (e,T;n)=e-^-*e(-T-- 



eo 



n \e F + e L 

The number of hot electrons N > . n (r) with energy greater than e F is 



N- 



n(r) = f 



f e {e,T; n)de = /o;„(ejr + e L ) 



1 + l 



-1 -l/n 



(10) 



(11) 



With j3 = X — 77+eT' ^ s m * nc ^ mear model, the number of hot electrons in the pile- up near the Fermi energy 

comes out to be: 



N P e;„(T) = / e 



r j 1 /n e -v(r+X n /n) 



J 7,0 



y " e v dy 



(12) 



The total number of hot electrons " Nh ot;n (T)" is iVft ot; „(r) = iV >;n (r) + N F;n (T). In the limit r << /?, or i << 
^ ( e ° ) . The expression for hot electrons in the pile up, and total number of hot electrons will reduce to: 



x -n _ v 



I _ e -(x -vh 



N h ot;n( T ) = /o( £ F + e L ) 



(1 + nx-V)- 1 /" + ^ (1 - e-fe" -'> T 



(13) 
(14) 



In the limit n going to zero, we recover the results of the linear model. 



III. INCORPORATING PUMP-PULSE DURATION 



Here, we consider the system being pumped by a rectangular femtosecond pulse of duration t p . The effect of the 
pulse can be taken as the convolution integral of the rectangular pulse with the respective hot-electron time-evolution 
curves(for the case n=0): 



1 r 

~* {T) = - L 

t p Jo 



i(t p ,t) 



Npe{r — x)dx, 



(15) 
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FIG. 3: Decay of hot electrons N pe - ;P (T) in the pile up as a function of r in the presence of a rectangular laser pulse with width 
t p = 3. Top most curve is for r\ — 0, lowest for 77 = 0.9 with a step of 0.1 




FIG. 4: Decay of hot electrons N hot - p (T) as a function of t, with same set of parameters as in Figure 3 



1 r-min{r p ,r) 

N p At) = - N hot (r-x)dx. 

T P JO 



(16) 



IV. CONCLUSION 



The distinctive feature of our calculated time evolution of the photoexcited non-equilibrium electron distribution is 
its peaking at the Fermi energy, as readily seen in Figure 1. It reflects the fermionic blocking effect. This can be, and 
should be, probed in a femtosecond pump-probe experiment. This calculation refers to a situation not describable 
by the usual two-temperature model. Here, instead, the electron system itself is comprised of two subsystems, one 
degenerate and the other nondegenerate. The nondegenerate subsystem relaxes towards the degenerate subsystem by 
energy loss to the lattice. 
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